On the scaling properties of 2d randomly stirred Navier— Stokes equation 
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nan's predictions [3, [5J, [21| . Very strong laboratory (see 
171 |26| and reference therein) and numerical evidences 



(see e.g. and references therein) support Kraichnan' 
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In two dimensions, the joint conservation of energy 
and enstrophy has relevant consequences for the Navier- 
Stokes equation. In 2d it is possible to prove that in 
the deterministic case the solution of the Cauchy prob- 
lem exists and is unique 12(1 and, very recently, that in 
the stochastic case H H M, El E| the solution is a 
Markov process exponentially mixing in time and ergodic 
with a unique invariant (steady state) measure even when 
the forcing acts only on two Fourier modes [13] • At large 
Reynolds number, the 2d scaling properties also differ 
from the 3<i-case. A long standing hypothesis [HI], very 
recently corroborated by numerical experiments [y, [jfl, 
surmises the existence of a conformal invariance in 2d. 
A phenomenological theory due to Kraichnan [l8| and 
Batchelor [2| predicts the presence of a double cascade 
mechanism governing the transfer of energy and enstro- 
phy in the limit of infinite inertial range. Accordingly, 
an inverse energy cascade with spectrum characterized 
by a scaling exponent dg = —5/3 appears for values of 
the wave- number k smaller than the typical scale Uf of 
the forcing scale. For wave- numbers larger than \zf a di- 
rect enstrophy cascade should occur. The corresponding 
energy spectrum scales as = — 3 + . . . where the dots 
here stand for possible logarithmic corrections. As em- 
phasized in 0, Q Kraichnan's theory is encoded in three 
hypothesis, (i) velocity correlations are smooth at finite 
viscosity and exist in the inviscid limit even at coinciding 
points, (ii) Galilean invariant functions and in particular 
structure functions reach a steady state and (iii) no dis- 
sipative anomalies occur for the energy cascade. Under 
these hypotheses, if the forcing field is homogeneous and 
isotropic Gaussian and time 5-correlated it is possible to 
derive asymptotic expressions of the three point struc- 
ture functions of the velocity field consistent with Kraich- 



s theory However, a first principle derivation of the sta- 
tistical properties of two dimensional turbulence is still 
missing. An att emp t in this direction has recently been 
undertaken 14, [lj| (see also 24 1 and [lij]) by applying 
a renormalization group improved perturbation theory 
(RG) ([27J and [l( for review of applications to fluid tur- 
bulence) to the randomly stirred Navier-Stokes equation 
with power law forcing. However how it will be discussed 
in details in the sequel, RG analysis leads to a scenario 
not a priori consistent with Kraichnan's theory. The goal 
of the present work is to shed light on this issue. 

RG starting point is the randomly stirred Navier- 
Stokes equation 

(d t + v- d)v a - v Q d 2 v a = -d a P + v a (1) 

d-v = 

where a — 1,2, P is the pressure enforcing incompress- 
ibility and £o is an Eckman type coupling providing for 
large scale dissipation. The forcing / is a Gaussian field 
with zero average and correlation 

-< / Q (x, t)f(y, s) y= S(t - s) F°"»( x - y) (2) 



Fq is a constant specifying the amplitude of forcing fluc- 
tuations and T a P(p) is the transversal projector. The 
function x in |4]) is slowly varying for m -C p -C M and 
set infra-red m and ultra-violet cut-offs M scales for the 
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forcing. Its detailed shape does not affect the scaling pre- 
dictions of the RG. Irrespectively of the spatial dimension 
d, the cumulative spectrum F°^(0) of the forcing (Ein- 
stein convection for repeated vector indexes) diverges for 
< e < 2 as the ultra-violet cut-off M tends to infinity 
hence providing for stirring at small spatial scales. For 
e > 2, F a a (0) is dominated by small wave-numbers and 
thus describes infra-red stirring. At e = 0, the scaling di- 
mensions of the material derivative, dissipation and forc- 
ing in JT]) coincide for a scaling dimension of the velocity 
field d v = 1 in momentum units. Thus e — provides a 
marginal limit around which scaling dimensions can be 
perturbatively determined with the help of ultra-violet 
RG. The main result 
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1 151 , Il6j is the existence of a 
non-Gaussian infra-red stable fixed point of the RG flow 
yielding for the energy spectrum the prediction 
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The adimensional function R depends upon infra-red 
scales m and kb = (e £o/-Po) 1 ^ 6 ~ 2e ' ) and admits a reg- 



ular expansion in powers of e at the RG fixed point 14 1 . 
The resummation leading to k 1 ^^ in ([5|) is derived from 
the solution of the Callan-Symanzik equation 27] which 
in the present case requires scaling at finite e to stem 
from the balance of the two terms in the material deriva- 
tive with the forcing i.e. from the requirement of Galilean 
invariance alone. According to [ll. ll4l . lla ] composite oper- 
ators at small e does not induce any self-similarity break- 
ing by the infra-red scales m and kb- The conclusion is 
that the spectrum should scale for small e with exponent 
rfg = 1 — 4e/3 as in the 3d-case [l|. Such conclusion 
seems at variance with Kraichnan's theory which instead 
suggests the occurrence of an inverse cascade at small 
e. Namely under the assumptions (i),(ii),(iii) [H, 0] the 
energy balance equation 

{dt +£o - wd 2 ) ■< v a ( x ,t)v a (0,t) y 

- l -d^< | \Sv\ | 2 (x, t) 5v» (x, t)y= F(x) (6) 

with£v a (x,i) = v a (x,t)-v a (0,t) yields for £ = 0, e < 2 
and r 2 := ||a;|| 2 



S 3 (r) ~ Cl F M 4 - 2e r + 



r 3-2s 



1 > r » M~ l (7) 



S3 denotes the three point velocity longitudinal structure 
function and Cj i = 1,2 two adimensional coefficients ir- 
relevant for the present argument. Power counting based 
on then predicts a dz = —5/3 inverse cascade spec- 
trum at small e with at most sub-leading corrections con- 
sistent with the RG prediction. 

The currently available numerical resources permit to 
contrast the two apparently discordant predictions ([3]) 
and ([7]) with the actual Navier-Stokes phenomenology. 



To attain this goal we integrated the Navier-Stokes equa- 
tion ([1]) for the vorticity field (u — € a pd^v a ) with a stan- 
dard, fully-dealiased pseudospectral method in a doubly 
periodic square domain of resolution 1024 2 . Time evo- 




FIG. 1: Energy spectrum for e = 0.5 and standard molecular 
dissipation. An inertial range with inverse energy cascade sets 
in for sufficiently small viscosity v = O(10~ 4 ) (solid line). At 
larger viscosity the inertial range is suppressed and dissipative 
spectrum £ ~ fc 1_2e is observed (dashed line). 



lution was computed by means of a standard second- 
order Runge-Kutta scheme. We repeated our numeri- 
cal experiments for different values of the hyperviscos- 
ity (— l) p+1 ^oi9 2p 7j including p = 1, standard viscosity 
as in ([T]), and p — 4, 6 obtaining qualitatively identical 
results. The outcomes evince that for all < e < 2 
and sufficiently small viscosity, an inverse energy cas- 
cade with exponent consistent with dg = —5/3 is ob- 
served (see Fig. [!}. If the viscosity increases, finite 
resolution effects prevent the observation of an inertial 
range as the Kolmogorov scale becomes of the order of 
the infra-red dissipation scale. In such a case a dis- 
sipative spectrum appears with scaling exponent con- 
sistent with the prediction dg = 1 — 2 e dictated by 
the balance between forcing and dissipation (see again 
Fig. [T]). In agreement with (J7J), local balance scaling 
becomes dominant in the range 2 < e < 3 (see Fig. 
[2]). There, energy spectra scale with an exponent com- 
patible with g?£ = 1 — 4e/3. Finally, for e > 3 a di- 
rect enstrophy cascade invades the whole inertial range. 
The phenomenology just described is confirmed by the 
inspection of the energy and enstrophy fluxes in the in- 
ertial range (see Fig. [3]). They are respectively defined 
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by U E (k) cx /* j^Vl -< ?{v} a (-«/)J{(« • 8v a )} (q) 

and U z (k) cc /* -< ?M (-g)J{(t> • du)} (q) >- 

where 5" denotes the Fourier transform and k a is cho- 
sen such that n^fc) and Tlz(k) are positive quantities 
in the inertial range. As spectra and structure functions, 
fluxes highlight the existence versus e of three distinct 
scaling regimes the origin of which can be understood 
by contrasting ^ with the energy and enstrophy inputs 
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FIG. 2: Third order structure function of longitudinal velocity 
increments S3(r) for e = 1 (squares), e = 2.5 (circles), e = 4 
(triangles). When e < 2 the presence of an inverse energy 
cascade with constant flux results in the scaling Ss(r) ~ r 
(solid line). When e > 3 the presence of a direct enstrophy 
cascade with constant flux results in the scaling 5*3 (r) ~ r 3 
(dashed line). For 2 < e < 3 the scaling is compatible with 
the RG argument (dash-dotted line) 

J B (fc) oc J ko dqqF v (q) and Jz(fc) oc J^dqqF v (q). For 
< e < 2 both Jg and Jz are peaked in the ultra-violet. 
Since stirring occurs mainly at small spatial scales an in- 
verse energy cascade sets in the whole inertial range. In 
the range 2 < e < 3 the energy input becomes infra-red 
divergent whilst Iz is still ultra-violet divergent in the 
absence of cut-offs. In this situation, the steady state is 
attained when the energy and enstrophy fluxes balance 
scale by scale the corresponding inputs (see Fig. [3]) . In 
this e-range the energy spectrum scales in agreement with 
the exponent which can be extrapolated but no longer 
fully justified using the RG. Finally, for e > 3 both the 
energy and enstrophy inputs become peaked in the infra- 
red. In such a case 0, 0] the right hand sides of Eq. 
© and of the analogous equation for the vorticity cor- 
relation admit a regular Taylor expansion for mx <C I . 
The hypotheses (i), (ii), (iii) thus recover S3 (x) oc r 3 i.e. 
Kraichnan's scaling for the direct enstrophy cascade in 
agreement with our numerical observations. 

As far as intermittency is concerned, our numerics 
support normal scaling of velocity structure functions 
(Fig. 3]). Higher order vorticity structure function are 
compatible with a weakly anomalous scaling (Fig. [5]). 
These results provide an indirect positive test for the oc- 
currence of a vorticity dissipative anomaly Q • 

In conclusion, our numerics fully support the validity 
of Kraichnan's theory for all values of the Holder expo- 
nent of power law forcing in 2d. RG scaling prediction 
is not observed in the range where it was supposed to 
appear. It seems to us that this in not trivially a conse- 
quence of low Reynolds number (Ofe 1 / 2 ) for Fq = 0(e) 
see discussion in section 9.6.4 of [131 ]) entailed by per- 
turbative expansion. Instead scaling seems to be related 




FIG. 3: The energy flux He (symbols) and the energy input 
Ie(Ic) (thin line) for different values of e. (a) e = 1; (b) 
e = 2.5; (c) e = 4. The energy input is defined as /^(fc) oc 
J * dq q F v (q) for e < 2 and as Ie (fc) oc dq q F v (q) for e > 2. 
In (b) the enstrophy input, Iz(k) oc f£ dqq 3 F(q) with and 
the enstrophy flux Ilz(k) are also plotted in order to show 
their balance at each inertial range scale. 



to the convolution of the response field with the forc- 
ing kernel in the Schwinger-Dyson equation 0, [2?| giving 
rise to a non-local scaling field of dimension do = —1/3 
independently of e. This point however deserves more 
theoretical inquiry. 
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FIG. 4: Structure function of longitudinal velocity incre- 
ments 5*4(r) (squares) and S§(r) (circles) vs. S2M. The lines 
represent the non- intermittent scalings 5*4 (r) ~ Siiff" (solid 
line) and S 6 (r) ~ S2(r) 3 (dashed line), (a) e = 1; (b) e = 2.5; 
(c) 6 = 4 
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FIG. 5: Structure function of vorticity increments Si?\r) for 
n = 1, 2, 3, 4 (square, circles, triangles, poligons). In the inset 
we show the scaling exponents £„. Here e — 4. 
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